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Possible Kondo effect in Pr skutterudite is studied with attention to characteristic features 
of low-lying crystalline electric field (CEF) levels and the conduction band. A mechanism 
for the small CEF splitting between a singlet and a triplet is proposed as combination of 
the point-charge interaction and hybridization of 4/ with ligand p states. Provided 4/ 3 con- 
figurations dominate over 4/ 1 as intermediate states, p-f hybridization favors the triplet, 
while point-charge interaction favors the singlet. For realistic parameters for hybridization 
as well as 4/ 1 and 4/ 3 levels, these singlet and triplet can form a nearly degenerate pseudo- 
quartet. It is found that one of two spin 1/2 objects composing the pseudo-quartet has a 
ferromagnetic exchange, while the other has an antiferromagnetic exchange with conduction 
electrons. The magnitude of each effective exchange depends strongly on a parameter char- 
acterizing the triplet wave function under the T/j symmetry. It is argued that differences of 
this parameter among Pr skutterdudites are responsible for the apparent diversity of their 
physical properties. Numerical renormalization group is used to derive the renormalization 
flows going toward singlet, doublet, triplet or quaret according to the CEF splitting and 
exchange interactions. 

KEYWORDS: crystalline electric field, hybridization, Kondo effect, filled skutterudite, 
PrOs 4 Sbi 2 , PrFe 4 Pi2 



1. Introduction 

Interactions between conduction electrons and localized degrees of freedom in metals can 
lead to profound effects due to multiple particle-hole excitations near the Fermi surface. The 
resistance minimum as a function of temperature was explained first by Kondo 1 who identified 
the crucial logarithmic term. In his model, a magnetic impurity is coupled antiferromagnet- 
ically with spins of conduction electrons. The essense of physics disclosed by Kondo 1 has 
far reaching consequences not only in condensed-matter physics, but also in particle physics; 
Kondo effect is the best example of the renormalization effect leading to the confinement and 
asymptotic freedom. 2 Since then many variants of Kondo effect have been found. We mention 
among others a case leading to a non-Fermi liquid ground state. 3 The main purpose of this 
paper is to show that a new aspect of Kondo-related physics appears in filled Pr skutterudite 
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compounds where astonishing varieties of anomalous properties have recently been discov- 
ered. The properties include heavy-fermion superconductivity and high-field ordered phase in 
PrOs4Sbi2, 4 ' 5 antiferro-quadrupole order in PrFe4Pi2, 6 metal-insulator and structural phase 
transitions in PrRu4Pi2- 7,8 While neutron scattering has observed clear CEF transitions in 
PrOs4Sbi2, 9 ~ 12 only broad quasi-elastic features are visible in PrFe4Pi2 above the tempera- 
ture of quadrupole order. 13 On the other hand, intriguing temperature dependence of CEF 
levels has been observed in PrRu4Pi2 below the metal-insulator transition. 14 For proper un- 
derstanding of these phenomena, interaction effects of 4/ electrons with conduction electrons 
should be taken into account. Especially challenging is to identify the parameters that control 
diversity of phenomena by their slight change. 

In this paper, we clarify the importance of hybridization between 4/ and ligand p electrons 
in giving rise to CEF splittings as well as possible Kondo effect. It is shown that the observed 
sequence of CEF levels can be understood as a result of competition between the point- 
charge mechanism and p-f hybridization. Furthermore, it is shown that one of two pseudo- 
spins forming the lowest four states has antiferromagnetic exchange with the conduction 
band. Kondo effect may appear strongly or only weakly depending on the wave functions 
and the level of the triplet CEF states. We argure that the delicate balance between the CEF 
splitting and Kondo effect should be the main reason for the apparently diverse properties of 
Pr skutterudites. 

2. Mechanism of CEF splitting in skutterudites 

2.1 Crystal structure and point- charge model 

In the filled skutterudite structure, each rare-earth ion is surrounded by 8 transition metal 
ions which form a cube, and 12 pnictogens which form an icosahedron deformed slightly from 
the regular one. The Pr site has local symmetry Th which has no four- fold rotation axis. 15 In 
this symmetry, 4/ 2 Hund's-rule ground states 3 H^ split into a Ti singlet, a T23 non-magnetic 
doublet, and two T4 triplets. Of these, T23 corresponds to T3 in the point group Oh, and two 
T4 triplets are written as and , which are linear combinations of T4 and T§ states in 
Oh- In addition to the Coulomb potential from ligands, hybridization between 4/ electrons 
and pnictogen p electrons contributes to CEF splittings. Under this crystal symmetry, the 
CEF potential is written as 



in the standard notation. 5 

In the point charge model, CEF coefficients A4 and Aq in eq.(l) are determined by coor- 
dination of the charge and the radial extension of the 4/ wave function. For transition ions, 



Vcef = MOl + 50 i] + A c 6 [O 6 - 210 4 ] + Al[Ol - 6 6 ] 



= W \x— + (1 - Ixl)-^- + v— 
x go + { i 1^1260 +y 30 



(1) 
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the eightfold coordination around Pr gives 

where f3j and 7j are Stevens coefficients, d t is the distance from the origin, and Z t is the 
effective charge per transition ion. There is no contribution to A\. On the other hand, the 
deformed icosahedron of pnictogens (X) gives 

A 4 = -2U(e )-^^(r*)(3j, (3) 
^ = -2/ 6 c (0 o )-^^<r 6 > 7 J, (4) 

where / n (0o) with n = 4, 6 are numerical factors dependent on the vertical angle 29q of the 
X-Pr-X triangle as follows: 

UWo) = \ (cos 4^o + cos 4 6 + sin 4 6 ), (6) 
/ 6 C (# ) = i(cos40 o + cos 4 9 + sin 4 6 - 11 sin 2 6 cos 2 0„), (7) 

/e(0o) = y (cos 6 - sin 6 O - cos 60 ). (8) 

To estimate the CEF potential, we use the lattice parameter in PrOs4Sbi2, 16 which gives d t = 
4.031, d p = 3.481 and #0 = 24.6°. By interpolating the data of ref.17, we use (r 4 ) = 3.4a| 
and (r 6 ) = 19ag with an the Bohr radius. With trivalent Pr, the effective charges satisfy the 
relation 

3 + 4Z t + 12Z P = 0. 

Figure 1 shows the CEF levels as a function of Z t . It is seen that with Z t > the singlet 
Ti is favored as the CEF ground state. For slightly negative Z t , the CEF potential becomes 
almost isotropic as pointed out by Takegahara and Harima. 18 The quasi-isotropy is due to 
high coordination numbers of pnictogens and transition ions, which tend to compensate the 
anisotropy. 

2. 2 CEF splitting by p-f hybridization 

Another important mechanism for CEF splittings is the covalency effect, or hybridiza- 
tion between localized and ligand orbitals. Before dealing with Pr skutterudites, let us first 
consider hybridization between 4/ n state and pnictogen molecular orbitals for general n. We 
follow the idea of Takahashi and Kasuya 19 who studied rare-earth monopnictides. The large 
Coulomb repulsion makes it reasonable to neglect the width of the relevant bands in the first 
approximation. Thus the band energy is replaced by its average. The hybridization is of the 
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form 

#hyb = 2 V T(a,p)Pl (a)ucT fr(f3)ua + ^C., (9) 
T-(a,l3)vo 

where the indices a, [3 distinguish different states with the same irreducible representation T. 
A member in the irreducible representation is specified by v, and a denotes a spin component. 

We treat hybridization by second order perturbation theory, and assume the following 
intermediate states: (i) 4/ n_1 and an extra electron in vacant states, and (ii) 4/ ra+1 and extra 
hole in filled states. Energy shifts in the case (i) are given by diagonalization of the following 
(2 J + l)-dimensional matrices: 

AE r (M',M) =J2^E- {a) (M',M), (10) 

a 

where M, M' are z-component of J, and w a denotes the weight of p component in the band 
with symmetry T(a), whose energy is represented by er( Q ). In the case (ii) we define the matrix 
A£tj (M f , M) in a similar fashion. Without multiplet splittings of intermediate states, we 
obtain 

ae-^m; m) = -a - „ rM) £E v t^Z1 

ua (3/3' A ™- 1 + £r («) 

x (M'\4 mu Jr((3)ua\M), (11) 
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A£+ (M', M) = -n r(a) V V V ^ V n^) 

x {M'\f np)va fl {p>a \M), (12) 

where A n ±i (> 0) are excitation energies to 4/ n±1 , and np( a ) is filling of the pnictogen T(a) 

state per spin. By using the commutation rule of / operators, we obtain the total energy shift 

by adding contributions from 4/™" 1 and 4/ n+1 intermediate states: 

— \V I 2 
AE r{a) (M', M) = -Y J . r( °'f' | 2n r(a) d r 8 MM , 
p A n+ i + |e r(a )| 



+ 



n T(a) 1 - «r(a) 



_A n+ i + |e r ( a )| A n _i + e r ( a ) 



x (^'l/J^^/roswIM), (13) 

where dr is the dimension of irreducible representation V. The first diagonal term does not 
contribute to the CEF splitting. 

For calculation of matrix elements in eq.(13) we need a matrix element of the type: 

(r-WL'M' L S'M' s \f mu(T \f n (LS)JMj}, (14) 

where we have rewritten \M) in more detail with M = Mj, and 7' is a label to distinguish 
states with the same L' and S' in 4/™ _1 . The matrix elements between states with plural 4/ 
electrons should respect antisymmetry of the wave function. This situation is taken care of by 
a factor (f n ^yLS{\f n ~ l j'L'S' , f), which is called the coefficient of fractional parentage. 21 We 
obtain 22 

(r-WL'M' L S'M' s \f mvo \r(LS)JMj) 

= (-ly-^r-YLSiir-WL's'j) 

x (L'M' L lm l \LM L )(S'M' s so\SM s )(lm l \lT{P)u)(LM L SM s \JMj). (15) 

miM L M s 

The matrix element in eq.(13) can be calculated with use of eq.(15) for any 4/ n_1 configura- 
tion. It is not necessary to compute the contribution of 4/ n+1 configurations separately. 

We apply this general framework to Pr ion in the point group O^. Figure 2 shows CEF level 
splittings of 4/ 2 states for different symmetries T2, T4 and T§ of hybridization. We have taken 
case of 4/ 3 intermediate states, and have assumed nr = 1/2 and ep = for each symmetry 
T. Since 4/ 1 and 4/ 3 intermediate states give opposite contributions to the level splitting, 
the sequence is reversed in the case of 4/ 1 intermediate states. In reality, the sequence is 
determined by competition between both intermediate states. 19 In eq.(13), the competition 
appears in the factor with np( a \. We have also done calculation where only Hund's-rule ground 
states in 4/ 2±1 are considered as intermediate states. The tendency of opposite contributions 
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Fig. 2. CEF level splittings from hybridization with 4/ 3 intermediate states in the point group Oh- 



from 4/ 1 and 4/ 3 states persist, although the relative distance is not completely the same. 



2.3 Slater-Koster parameterization 

We proceed to the actual Th point group around Pr ion in which 4f electrons hybridize 
with pnictogen molecular orbitals with the symmetry a u (= T2) and t u . Hereafter we adopt 
the Mulliken notation such as a u for orbital symmetry, and the Bethe notation such as Pi for 
a representation with spin-orbit coupling. Figure 3 shows relevant molecular orbitals schemat- 
ically as a combination of atomic p orbitals at pnictogen sites. The hybridization is parame- 
terized by the Slater-Koster parameters (pfcr) and (p/vr) 23 as shown in Table I. It turns out 
hybridization with a u comes only from (pfir). 

According to band calculation, 24 the conduction band is formed mainly by a u orbitals. 
On the other hand, two t u bands are a few eV above the Fermi level and three are a few eV 
below. We neglect contributions from t u bands to CEF splittings for simplicity. We expect the 
following argument still applies to real Pr skutterudites at least qualitatively by the following 
two reasons: Firstly, the effect of t u bands is small as compared with the a u band because 
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Table I. Hybridization constant Vr( a ,p) between / and p states with the same symmetry. The angle 
6*o has been defined below eq.(5), and s = sin9 ,c = cos^o- 

of the larger excitation energy. Secondly, since 4/ 1 and 4/ 3 intermediate states give opposite 
contributions to level splittings, the effect of two empty and three filled bands tends to cancel 
each other. 

With only a u band taken into account, eq.(13) shows that the CEF level sequence is 
determined by two parameters; occupation of the a u band and A3/A1. If the latter is less 
than unity, the sequence becomes the same as that in 4/ 3 intermediate states. Otherwise 
the sequence is reversed. To a good approximation we can assume the half-filled a u band, 
and er = for the band. According to available information on PrP, 19 we obtain the ratio 
A3/A1 = 0.6. Hence the level structure due to hybridization alone is dominated by 4/ 3 
intermediate states, and should be given as shown in the leftmost part in Fig. 2. However, 
we have to consider both point-charge interaction and hybridization for the final sequence of 
CEF levels. 

2.4 Combination of Coulomb interaction and hybridization 

We combine both contributions to CEF splittings by point-charge interaction and by 
hybridization. Figure 4 shows computed results as a function of strength of hybridization, 
(p/7r) 2 /A_ where we have introduced 

1/A_ = I/A3 - 1/Al 

For the point charge parameters we tentatively take Z t = 2 and Z p = —11/12. We note 
that in Mossbauer experiment on PrFe4Pi2, 20 Fe is reported to be trivalent. On the other 
hand, the energy-band calculation suggests a smaller valency for Fe. 24 Overall CEF splittings 
become larger with larger Z t as seen from Fig.l. As hybridization effect increases, r^ 2 -* level 
is stabilized and crosses the Ti ground state. The parameter (p/7r) 2 /A_ = 190K is consistent 
with the combination of Ai ~ 5eV, A3 ~ 3eV and (pfir) ~ 0.35eV. The (pfir) is about 50% 
larger than that obtained for Pr monopnictides. 19 Since the distance d p is smaller than the 
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Pr-P distance of 4.2A, the larger magnitude is reasonable. The level repulsion between 
and rjj 2 ^ around (p/7r) 2 /A„ = 110K is due to mixing of wave functions, and is a characteristic 
feature in the point group T^. In contrast, a level crossing occurs in the Oh group since there 
is no mixing between T4 and T$. 

3. Exchange interaction in the pseudo-quartet 

We derive the low-energy effective interaction between 4/ and conduction states by re- 
stricting to the singlet and the lowest triplet for the CEF states. We shall show that the mag- 
nitude of the exchange interaction depends strongly on the nature of triplet wave functions 
even though the strength of hybridization and the CEF splitting are similar. This difference 
of effective exchange can be the main reason for the diversity of physical properties in Pr 
skutterudites. Although exchange interaction for f 2 systems has already been discussed in 
literature, 25 ' 26 the feature mentioned above has so far escaped theoretical attention. 



3.1 Second- order exchange interaction 

For simplicity, we neglect the multiplet splittings in intermediate states with 4/ 1 and 4/ 3 
configurations. Only the a u state is kept for conduction band, which intersects the Fermi 
surface and is nearly half-filled. We take the hybridization parameter V-2 U real and deal with 

H hyh = V2 U J2f& + h.c, (16) 

<j 

where c a annihilates a conduction electron in the Wannier orbital with a u symmetry at the 
origin, and /«J creates an 4/ electron with the same (a u ) orbital symmetry. In the second-order 
perturbation theory, the effective interaction is given by 

#int = PH hyh — —QH\ lyh P, (17) 

where P is the projection operator onto 4/ 2 states, and Q = 1 — P. We have already taken 
into account such part of the second-order hybridization that is diagonal with respect to the 
CEF states. Now we deal with off-diagonal components. 

Let us fist consider 4/ 1 as intermediate states, and neglect the multiplet splittings. Then 
we obtain 

1 MM' acr' 

where M denotes azimuthal quantum number of J = 4, and we have introduced the quantity 
A(MM';aa')= £ {M^rmm^immMM'). (19) 

mira s 

Here mi and m s denote azimuthal quantum numbers of / = 3 and s = 1/2, respectively. 

In the case of 4/ 3 intermediate states without multiplet splittings, the effective Hamilto- 
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nian is given by 

#i„t[4f 3 ] = E E A'(MM'-aa') \M)(M'\ C t C(T ,, (20) 

3 MM' crcr' 

where we have introduced the quantity 

A'(MM';aa') = E(M|/ CT |/ 3 )(/ 3 |/^|M'), (21) 
/ 3 

with summation over all f 3 states. Because of completeness of the intermediate states, we 
obtain the relation 

A(MM'; aa') = S M M'S a a' ~ A' (MM'; aa'), (22) 

where the first diagonal term is irrelevant. 

Hence consideration of 4/ 3 intermediate states in addition to 4/ 1 is accomplished by the 
replacement 

V 2 2 JA l - VKl/A, + 1/A 3 ) = Vl/A, 

in eq.(18). Note that 4/ 1 and 4/ 3 states contribute additively in contrast with the case of 
CEF splittings. Although the bare diagonal terms are already included as CEF splittings, 
higher order potential scattering terms contribute to renormalization of CEF levels as well as 
damping of states. 

3.2 Symmetry analysis in the pseudo-quartet 

Let us first write down the CEF eigenfunctions in Oh in terms of J = 4 states as 

|ri) = ^l+4) + ^|0) + yi|-4) > (23) 

|r 4 ,±) = T yi| T 3) T y / ||±i), (24) 

|r 4 ,0) = y|(|+4)-|-4)), (25) 

|r5 ' ±) = ± \/I l±3)T Vf lTl) ' (26) 

|r 5 ,0) = y|(|+2)-|-2)). (27) 

Following ref.27, we represent the lower triplet eigenfunctions in Th in terms of the functions 
obtained above: 

\T t ,m) = y/w\T^m) + y/1 - u;|r 5 ,m), (28) 

where the index m = 0,±1 distinguishes the expectation value of J z , and the parameter 
< w < 1 measures the weight of ^4). The other linear combination with a negative coefficient 
has a higher energy in our CEF potential. 
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(29) 



where r 4 has the same symmetry as the magnetic moment. Thus we see that both CEF states 
are magnetic. On the other hand, the direct products of pseudo-quartets can be decomposed 
as 



(Ti © r 4 ) © (ri © r 4 ) 
(r 1 ©r 5 )©(r 1 ©r 5 ) 



ri + 2r 4 + (ri + r 3 + r 4 + r 5 ) 
ri + 2r 5 + (ri + r 3 + r 4 + r 5 ) 



(30) 
(31) 



Hence, the cross product of Ti and T 4 produces newly two T 4 representations, of which one 
is time-reversal odd and the other is even. The odd representation represents the magnetic 
moment, while the even one the hexadecapole moment. On the other hand, the pseudo-quartet 
Ti © ]?5 does not produce a new T 4 representation. Physically, this means that T 4 as the first 
excited states gives rise to a van Vleck term in the magnetic susceptibility, while Ts does not. 

It is convenient to introduce the effective angular momentum operators in the pseudo- 
quartet in T/j. In the case of w = (pure Ts triplet), the magnetic moment within is the 
only relevant quantity. The vector operator in this case is written as X t . On the other hand, 
another vector operator X s is necessary to describe the magnetic moment of van Vleck type. 
In the matrix representation, we obtain 
/ 
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where the basis set is arranged in the order of 
rules of these operators are given by 



|r t ,+),|r t ,o),|r t ,- 



[ X H X j] = i e ijkX k , 

[Xf,Xj] = ieijkXl, 
[Xj,Xj] = ieijkXl, 



TrXjX] 



2,0<yy/5ij , 



( 1 " 



10 
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(33)' 
The commutation 

(34) 
(35) 
(36) 
(37) 
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where 7, 7' represent either t or s. 

We project the basis set \M) to the pseudo-quartet. In terms of two vector operators X t 
and X s , we obtain the following form for the effective interaction within the pseudo-quartet 

#int = ^ E ( atXt + asXS ) • ^cic^, (38) 

a/3 



where 



10 + 88u; 4 5w , on , 

at = -^155-' as = -33VT- (39) 
In deriving eq.(39), we used the relation 

{M\fl,\mim s ) = -y/2 Yl {JM\LM L SM S ) 

m[M L Ms 

x (LMLllmilm'^iSMs \ sm s sa'} {lm\ \IT 2 ), (40) 

with (lmi\lT2) = ±l/y/2 (for mi = ±2). Both a t and a s are negative, which means ferromag- 
netic exchange between conduction and 4/ moments. The origin of the ferromagnetic sign 
for a t is traced to the Hund rule involving the spin-orbit interaction; the dominant orbital 
moment is pointing oppositely from the spin moment. In other words, the spin exchange is 
antiferromagnetic as is usual for a hybridization induced exchange. The sign of a s , on the 
other hand, is not physially meaningful. 

The relation to the pseudo-spin representation of ref.27 is as follows: 

X 1 = Si + S2, X s = Si — S2, 

where Si and S2 are spin 1/2 operators. It is most natural to interpret the commutation rules 
eq.(37) in terms of Si and S^. The exchange interaction in terms of pseudo-spins is given by 

H ink = (J 1 S 1 + J 2 S 2 )-8 c , (41) 

where Ji = (a t + a s )V 2 2 u /A, J 2 = (a t - a s )V 2 2 u /A and 



u 2 

a/3 

Figure 5 shows the coefficients a t ±a s of pseudo-spins. It should be noticed that a t — a s becomes 
positive for w > 0.00324, and represents antiferromagnetic exchange J2. The emergence of 
antiferromagnetic exchange is due to the particular CEF level structure in Pr skutterudites. 
The antiferromagnetic exchange interaction measured by a t — a s is almost negligibe in the 
pure T5 case (w = 0), and becomes an order of magnitude larger as w increases toward unitiy, 
i.e., toward pure It is likely that the Kondo-type behavior seen in PrFe4Pi2 originates 
from a large value of w together with a small singlet-triplet splitting. 
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4. Renormalization group analysis of the singlet-triplet Kondo model 

In this section we classify all fixed points of the model characterized by the exchange 
interaction of eq.(41). Although signs of J\ and J2 are restricted in Pr skutterudites, we 
consider here all combinations of signs. Possible Kondo effect in systems with a CEF singlet 
ground state has been discussed from various aspects in literature. 28-30 Our model relevant 
to Pr skutterudites has a unique feature of single channel of conduction electrons plus two 
localized spins with S = 1/2. It will be shown below that the single-triplet CEF system has four 
fixed points with different residual entropies. These fixed points are understood naturally only 
if we consider renormalization of CEF levels. Hence we add the CEF terms to the interaction 
given by eq.(38) and obtain 

H s . t = e t Pt + e s P s + (/ t X* + I S X S ) • s c , (42) 

where It = atV^/A, I s = agV^/A and where et and e s describe the CEF levels. Only the 
difference Acef = ft — £s is relevant. With use of eq.(41), the interaction is alternatively 
written as 

H s . t = A CEF Si • S 2 + (J1S1 + J2S2) • s c . (43) 

We call this model, combined with kinetic energy of conduction electrons, the singlet-triplet 
Kondo model. The model has interesting features of renormalization because the screening 
channel is single (a u ), while localized spins are two kinds. 
The obvious fixed points of the model are the followings: 

(i) CEF singlet in the case of Acef > | Ji \,\ J2 \ , 

(ii) CEF triplet in the case of -A C ef > l^l, 

(iii) doublet ground state in the case of — J\ J2 S> | Acef | , 

(iv) quartet ground state in the case of both J\ and J 2 being renormalized off, and Acef = 0. 
The fixed point (iv) can be regarded as a special case of (i) or (ii) and is accidental. The nature 
of the fixed point (iii), where both the ferromagnetic J\ and antiferromagnetic J2 exchanges 
are present, can be understood as follows: The component S2 will make a singlet by Kondo 
screening, while the ferromagnetic exchange of other component S± will be renormalized off. 
As a result, the fixed point will be the doublet made up of S\. 

We have performed explicit calculation using the numerical renormalization group (NRG) 
for this model. Figure 6 shows the entropy against temperature in the case of Acef = for 
various cases of exchange interactions. We used the starting parameters J\/D = ±0.8,0 and 
J2/D = ±0.44,0 where D is the half- width of the conduction band with constant density 
of states. In Fig.6 the entropy in the case of J1J2 < (AF,F) goes to zero. Since we take 
Acef = 0, one might naively expect the fixed point (iii) in this case with the residual entropy 
In 2. The computed result is naturally interpreted as renormalization of Acef to a positive 
value. Namely the ground state is connected with the CEF singlet. 
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In the case of J\ < 0, J2 < (F,F), the fixed point becomes the triplet (ii), instead of 
quartet (iv) as might be expected naively. Again the result is naturally interpreted as renor- 
malization of Acef which now flows to a negative value. If one started with the interaction 
term such as eq.(41) without Acef, the renormalization flow to triplet or singlet would seem 
mysterious. 

There is a tiny critical interaction J2 < at which the singlet fixed point is taken over 
by the doublet fixed point. The critical interaction gives an unstable fixed point at which 
non-Fermi liquid behavior should be realized. In contrast with the two-impurity Anderson 
model, we do not have a state connecting to two Kondo singlets since one of the exchange is 
ferromagnetic. Further details of the NRG results will be discussed elsewhere. 

In our model of Pr skutterudites with J\ J2 < (AF,F), the renormalization favors the 
singlet ground state. Even though the triplet CEF level is lower in second-order perturbation 
theory, higher-order terms may reverse the levels. If the system undergoes the metal-insulator 
transition, as in PrRu4Pi2, the singlet may lose the relative stability because the energy gap 
created in the conduction band will reduce the renormalization effect. As a consequence, the 
triplet CEF can be more stabilized as the energy gap increases. This situation may be relevant 
to the apparent reversal of CEF levels reported in PrRu4Pi2- 14 In the next section, we shall 
discuss another possibility for the temperature-dependent change of CEF levels. 

5. Discussion 

5. 1 Comment on other mechanisms for CEF splittings 

In addition to the point-charge mechanism and hybridization, there are other sources for 
CEF splittings such as anisotropic Coulomb and exchange interactions. 19 In the case of rare- 
earth monopnictides, these other effects involve 5d electrons and tend to cancel each other. 
In Pr skutterudites, there is no reason to expect enhancement of anisotropic Coulomb and 
exchange effects, since the conduction band is composed mainly of p electrons of pnictogens. 
Thus we expect that our analysis with two main effects shows realistic origins of CEF splittings 
in skutterudites. However, we have totally neglected intersite interactions of 4/ electrons which 
are essential to ordering phenomena. This aspect should be studied on the basis of CEF level 
structures obtained here. 

5. 2 Interpretation of different behaviors in actual Pr skutterudites 

It is known that the magnetic easy axis is different in PrFe4Pi2 and PrOs4Sbi2- Thus the 
CEF sequence seems also different in both cases. In the paramagnetic phase, PrFe4Pi2 has 
a sequence of magnetization M(100) > M(110) > M(lll) along each direction of magnetic 
field. 6 The anisotropy in PrFe4Pi2 is consistent with CEF transitions between T4 and T% in 
the notation of Oh group. Note that the van Vleck term from Ti — T4 transition is isotropic. 
Then candidates for the CEF ground state are Ti and T3. 6 To explain the first-order transition 
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to the quadrupole order, Kiss and Fazekas proposed a model where T\ is the ground state. 31 
The sequence of CEF levels appears like the left end of Fig. 4 where T3 is written as in the 
notation of T^. In the following discussion we assume the Kiss-Fazekas model for PrFe4Pi2- 

On the other hand, PrOs4Sbi2 shows the magnetic anisotropy M(110) > M(lll) > 
M(100) in the paramagnetic phase. 32 The anisotropy is consistent with CEF transitions be- 
tween and T3 in the notation of Oh group. The sequence near (p/7r) 2 /A_ ~ 180K in Fig. 4 
is consistent with this magnetic anisotropy. In fact, neutron scattering observed transitions 

from not only from the T\ level but from the first excited triplet T^. 12 The overall intensity 

(2) 

suggests that the weight w is small in this T\ states. 

By comparing these facts in PrFe4Pi2 and PrOs4Sbi2, we conclude that CEF parameters 
as shown in Fig. 4 vary substantially in Pr skutterudites. Experimentally, hybridization effect 
appears stronger in PrFe4Pi2 than in PrOs4Sbi2- At first sight, it is paradoxical that CEF 
levels in PrFe4Pi2 look similar to the one in the left end of Fig. 4. We suggest that the variety 
of CEF parameters comes from strong cancellation effects between hybridization with 4/ 3 
and 4/ 1 intermediate states. If one assumes nearly equal contribution from both intermediate 
states in PrFe4Pi2, we obtain small |1/A_| = |l/Ai — I/A3I. Thus even though V 2 m itself 
is large, the CEF levels may be similar to that given by the point-charge model. As 4/ 1 
intermediate states become more dominant than 4/ 3 states, T3 level becomes stable. This 
situation corresponds to negative Ai, or reversal of levels in Fig. 2. If such is a case, PrFe4Pi2 
may have T3 as the CEF ground state. 

Available information on the ratio A3/A1 in actual Pr systems is as follows: Pr 4/ spec- 
trum has been observed in X-ray photoemission spectroscopy (XPS) of PrFe4Pi2- 33 ' 34 How- 
ever, bremsstrahlung isochromat spectroscopy (BIS) for Pr skutterudites has not been carried 
out to our knowledge. Data are available for rare-earth metals such that A n+ i in Ce, Pr and 
Nd are 3.5, 2.1 and 1.7 eV, respectively, and A n _i are 0.3 ~ 1.9, 3.3 and 4.7 eV in the same se- 
quence. 35 In the case of monophosphides, it is reported that A n+ i/A n _i of CeP, PrP and NdP 
are 2.5, 0.6 and 0.5, respectively. 19 The tendency that A n+ i/A n _i decreases as n increases is 
also seen in light rare-earth hexaborides. 36 Then it seems reasonable to assume A3/A1 ~ 0.6 
in Pr skutterudites. Of course, hybridization with t u ligand orbitals may contribute to the 
competition to some extent although they were neglected in the present analysis. 

Unusual temperature dependence of CEF levels has been found by neutron scatering in 
PrRu4Pi2 below the metal-insulator transition temperature. 14 In the insulating phase, there 
appear two inequivalent Pr sites Prl and Pr2 with different CEF splittings. The observation 
of ref.14 can be qualitatively accounted for if we assume that (p/7r) 2 /A_ increases as temper- 
ature decreases below the transition point. To illustrate the idea, let us suppose in Fig.4 that 
(p/7r) 2 /A_ ~ 100K above the metal-insulator transition. Then suppose that (p/7r) 2 /A_ in- 
creases by about 10 K in Prl sites, and by about 100 K in Pr2 sites. As a result, Pr2 sites have 
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the triplet CEF ground state, while Prl sites keep to the singlet. However, the triplet changes 
the character substantially by the small shift of (p/7r) 2 /A_. It is an interesting open problem 
to identify the basic mechanism that causes the change of (p/7r) 2 /A_ in the insulating phase. 

5. 3 Dependence of Kondo effect on CEF wave functions 

Even with the same magnitude of Viu, the strength of exchange interactions can be much 
different depending on the CEF wave functions. As shown in Fig.5, the exchange interactions 
become weak for small w, namely for states close to r§ states. This is consistent with the fact 
that PrOs4Sbi2 does not show Kondo effect in resistivity. 4 On the other hand, if the CEF 
triplet consists of states with w ~ 1, i.e., close to T4 in Oh group, the exchange interactions 
can be strong. This is again consistent with the CEF states suggested for PrFe4Pi2, which 
shows Kondo effect in resistivity, 37 and where CEF excitations are hardly visible by neutron 
scattering in the paramagnetic phase. 13 It should be interesting to derive resisitivity and 
dynamical susceptibility theoretically for different values of w. The results of such study will 
be reported in the near future. 

6. Summary 

We have investigated the CEF level structures in Pr skutterudites with particular attention 
to the fact that hybridization through 4/ 1 and 4/ 3 intermediate states give opposite contri- 
butions to CEF level splittings. Provided 4/ 3 intermediate states are dominant, hybridization 

(2) 

with the a u conduction band stabilizes T\ triplet. On the other hand, positive point charges 
on transition metals stabilize T\ singlet. Thus a competition between point charge potential 
and p-f hybridization can lead to the experimentally proposed level structure for PrOs4Sbi2, 

(2) 

i.e., Ti ground state and T 4 low lying excited states. The difference between the CEF level 
structure in PrOs4Sbi2 and that in PrFe4Pi2 is ascribed to difference in competition between 
4/ 1 and 4/ 3 intermediate states rather than difference in strength of hybridization. It is shown 
that importance of Kondo effect depends strongly on the character of the CEF triplet next 
to the singlet. The predicted tendency according to each CEF wave function is consistent 
with absence of Kondo effect in PrOs4Sbi2 in resistivity and strong logarithmic anomaly in 
PrFe 4 Pi2- 
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Fig. 3. Illustration of pnictogcn molecular orbitals. Light and dark parts show the sign of p-wave 
functions at pnictogen sites. 
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Fig. 4. CEF parameters W,x,y and level structures derived from hybridization and point charge 
potential as a function of (p/7r) 2 /A_ . Point charges on transition metal and pnictogen are assumed 
to be +2 and —11/12, respectively. 
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